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We present a theory of Landau-Lifshitz-Gilbert damping a for a localized spin S in the junction coupled to the 
conduction electrons in both leads under an applied volatege V. We find the voltage dependence of the damping 
term reflecting the energy dependence of the density of states. We find the effect is linear in the voltage and 
cotrolled by particle-hole asymmetry of the leads. 

PACS numbers: 75.80.+q, 71.70.Ej, 77.80.-e 



INTRODUCTION 

Spintronics is an emerging subfield that holds the poten- 
tial to replace conventional electronic devices with spintronic 
analogues where the manipulation, control, and readout of 
spins will enable novel functionality with no or little elec- 
tronic charge dynamics In order to realize this promise, 
the spin dynamics of the small scale devices needs to be well 
controlled. One of the most pressing questions concerns a set 
up which would preserve coherence and allow a manipula- 
tion of spins. In most systems, the relevant spin degrees of 
freedom are coupled to some bath, such as a fermionic bath 
of electrons. The detailed dynamics of single spins when in 
contact with such a bath plays a pivotal role in addressing de- 
coherence in potential spintronic systems. 

The conventional way to treat this problem is via a 
Caldeira-Leggett approach where the external bath is modeled 
by collective excitations which are capable of destroying co- 
herent spin dynamics. Often, spin dynamics is described by a 
Landau-Lifshitz-Gilbert equation |2j,|3[: 
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both correlations. In the current work, we examine the volt- 
age dependence of Gilbert damping. For the sake of clarity, 
we depart from the Keldysh contour formalism of |71 |gl fcj], 
and use a Caldeira-Leggett approach. 

In what follows, we consider the case of a junction between 
two electrodes that contains one spin 5, see PigEJ This spin 
5 may be the spin of a single magnetic impurity or it may 
portray the spin of a cluster at low temperature when the spins 
in the cluster are locked. Upon applying a finite bias between 
the electrodes of FigfJ a current flow is generated. Thereafter, 
at long times, the system is at a steady but non-equilibrium 
state so long as the voltage bias V is applied. We will focus on 
the voltage dependence of the damping term a(V) in Eq.(l). 
We find that the change in the density of states associated with 
the chemical potential gradient across the junction triggers a 
modification to the damping a that is linear in voltage and is 
proportional to the particle-hole assymmetry of the density of 
states. The scale of the correction is set by the Fermi energy 
of the metal in the leads Ep and by particle-hole asymmetry 
in the density of states: 



(1) 



a(V)=a + a 1 (V) = a (l + O{l)eV/E F ). 



(2) 



where h is, up to constant prefactors, the external magnetic 
field and the coefficient a captures the damping due to the 
external bath. A caricature of the solution of this equation 0] 
is provided in Fig^ There are standard methods to calculate 
a in an equilibrium situation when, say, one considers a spin 
in a Fermi liquid 

In the current publication, we address a related novel ques- 
tion concerning the effect of an applied voltage bias on the 
Gilbert coefficient a. Our work complements the recent re- 
sults of 0] wherein the effects of the "retarded" electronic 
contributions in the equations of motion for a system of spins 
were studied. Both such retarded correlations |8| as well as 
additional "Keldysh" correlations generally manifest them- 
selves in the single spin equations of motion, see e.g. 0] 
for general spin equations of motion entailing the effects of 



This result vividly illustrates the presence of voltage induced 
damping in such junctions. Spin unpolarized electrons tunnel- 
ing across the junction interact via exchange interaction with 
the spin 5 and produce random magnetic fields that disorder 
the local spin. This noise augments that already present equi- 
librium magnetic noise in a Fermi liquid bath. Such a behavior 
of Oc(V) with the external voltage is in line with the works of 
flfjll . An analysis of a related single spin problem in a Joseph- 
son junction (instead of the normal junction studied here) was 
advanced in (111 EH. 

We will shortly derive the effective single spin action from 
which the principle equation of motion of Eq.lTiJ follows. Sev- 
eral technical details of our derivation are given in the appen- 
dices. 
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Hamiltonian governing this system is given by 




(a)a=0 



(b)a>0 



FIG. 1: Sketch of the dissipative spin dynamics. Panel (a) depicts 
a cartoon of the Larmor precession of the spin about the direction 
of an applied magnetic field (B). In panel (b), a caricature of the 
spin dynamics in the presence of Landau-Lifshitz-Gilbert damping 
is shown. 




FIG. 2: Magnetic impurity coupled to two electrodes, [j^ and /jr de- 
notes the chemical potentials of the left and right leads respectively. 
The voltage drop across the two electrodes is eV = MR- 



THE SYSTEM AND OUR PRINCIPLE RESULTS 

The physical system under consideration in this publication 
is illustrated in Fig|2] It consists of two (left(L) and right(R)) 
electrodes across which a voltage bias is applied; a magnetic 
impurity (S) is situated in between (or lies on one of) the elec- 
trodes. An external magnetic field B is present. In the absence 
of effects stemming from conduction electrons in the tunnel- 
ing barrier, the single spin would precess at the Larmor pre- 
cession frequency about the applied field direction (see, e.g. 
panel (a) of Fig^. With the external circuit elements present, 
the spin motion becomes dissipative (as schematically shown 
in panel (b) of Fig^. 

With the spin embedded in the tunneling barrier, the work 
function is modified and the conduction electron tunneling 
matrix element is supplanted by a Kondo like exchange term 
J {§ -Oc), with a c denoting the conduction electron spin. In 
what follows, we will dispense with the c subscript. The 



9t 
A 

A 
Ar 



?{ e + !H: s + !f{T, 

52 ^0-k c aka Calca > 
a,k,a 

-h-S(t), 



(3) 



a 



(4) 



a,&,op.p,<y 



where cL.„ (c a/ t ) creates (annihilates) an electron with mo- 
mentum k and spin a e {1,1} m the lead a 6 {L,R}. The 
abbreviation = e a ^ — /J a , where e a/ t is the energy of elec- 
tron with momentum k in (the lead) a and /j a is the chemical 
potential in (the lead) a. The second term in Eq.(@J, H s , is 
Zeeman energy of the spin in an external magnetic field B. 
Here, h = g^gB with g gyromagnetic ratio and /jg the Bohr 
magneton. The last term in Eq.l|4j, ttj, represents both Kondo 
coupling and direct tunneling process, where the amplitudes 
{ T ap} = {TllAlrArlArr} are tunneling matrix elements 
and their explicit forms are 



Tll = Ju.@-S(t)), 
T RR = J m (a-S(t)), 
Tlr = Trl = (7b8 ao / 



-■MS-S(O)). 



(5) 



Here, Tq is the direct tunneling matrix element and / a p is the 
Kondo coupling, 'while Q. denotes the Volume of each lead 
(assumed, for simplicity, to be the same)'. Typically, from 
the expansion of the work function for tunneling, Jlr/Tq ~ 
J/U, where U is the height of a spin-independent tunnel- 
ing barrier and / the magnitude of the spin exchange inter- 
action fi3l . Typical values of the ratio the spin dependent 
to spin independent tunneling amplitudes in Eq.(|5}, Ap/Ao, 
are 0(10 ), with a typical Fermi energy E§ of the order of 
several electron-volts. From the Hermiticity of the Hamilto- 
nian, we can find that the matrix element (TIxp)oo' satisfies 
((TapW)* = (Tpa)o'a- 

In the up and coming, we derive the effective action for 
the single impurity spin via an imaginary time path integral 
formalism. The full action is given by 



^ E 4to 3 x c «to + iS(d(S(x)) + \ dxtt (t), (6) 



where the second, Wess-Zumino-Novikov-Witten(WZNW), 
term in Eq.(|6j depicts the Berry phase accumulated by the 
spin fl3ll . In our action, we have the following quadratic form 
of fermions, 

p i 

dZ a £ £ c ato( 8 aP 8 °o'( 3 x + ^ak) + AapW)cppa' 
" atoppo' 

/•P 1 

= / rfX oEE W( M «p( T ) W)*P %><?'• (7) 
70 "atoppo' 
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We may integrate over the lead electrons to obtain the effective 
action for the spin 



Gilbert coefficient a in Eq.Q is exactly equal to r|. By vary- 
ing the total action with respect to the spin 8S/8S(x) = 0, we 
immediately obtain the Landau-Lifshitz-Gilbert equation with 



Seff(SM) ~ iS(i>(S(z)) + / dzh ■ S(z) - IndetM, (8) a = f A»)(see APPENDIX C). In other words, the voltage 
Jo dependence of a in Eq. (fill is identically the same as that of 



where detM means functional determinant of M. From the 
third term in Eq.©, we obtain a quadratic non-local in time 
interaction of the spin with itself, S(x) as 



A5(S(x)) = -2 dz dz'S(z) ■S(z')K(z-z') 



(9) 



where 



K(z) = £ * a p(x), 
a,pe{L,R} 



K, 
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with /(£) denotes the Fermi distribution function (see AP- 
PENDIX A). The effective action AS of Eq.© can be de- 
composed into two (trivial and non-trivial) components as 
AS = ASi c + ASdis, with 

ASi oc = -2K{(Q = 0) J dz(S{z)) 2 , 

AS dis = I dz [ dz'{S(z)~S{z)) 2 K(x-z'). (11) 



Here, K{ia = 0) is the zero-frequency Fourier component of 
K(z). The first term (A5i oc ) is a trivial constant as S(z) 2 = S 2 . 
The nonlocal part (A5dis) represents the dissipative effect due 
to the coupling between S(x) and electrons bath. The integral 
kernel K(z) is calculated in the same way as the Caldeira- 
Leggett theory OEl leading to 

g( x)=r^( a ) gg^#g > (i2) 

w Jo sinh(pco/2) 
where /(co) is the spectral density and its explicit form is 
/ap/po 



J ( & ) = E E E ^ - /M«(<» + - ^ 

apfceapep 



(13) 

The details of the derivation of Ea. (ll3> are provided in AP- 
PENDIX B. The spectral density of Eq.flSl. /(co), is esti- 
mated as 

/(a) = £ 4^ r^ &) r^pJV(^) 



[/(^)-/(^p)]5(fi)-^-^) 



■AxpJpc 



ap 



iV(^„ = 0)^=0)0), (14) 



where £" denotes the Fermi Energy of the lead a. It is obvious 
that /(co) in Eq. (I14> is proportional to co, i.e., /(co) is Ohmic. 
If spectral density is expressed as /(co) = T|Co/27C, then the 



dependence of a in Eq.Q is identically the same as that of 
/(co). We next examine the voltage dependence of /(co). 

If we apply a voltage leading to a chemical drop of Afj F — 
A/j r = eV. Assuming, for example,that the net charge on both 
right and left leads is unchanged, we also have D L {Ef )Afii + 
D R (Ef )A/jr = 0. With these constraints we get 



Aul = 

A/JR = 



D R (E F ) 



D L {E F )+D R {E F ) 

•f) 



D L (E F 



D L (E F )+D R (E F ) 



eV, 
eV, 



(15) 



the Gilbert coefficient a may be approximated as 



a(V) ~ 2ni^fD L (E F + Ap L ) 2 + J -fD R (E F +Ap R ) 
+jl R D L (E F +A r u L )D R (E F +A r j R ) 
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[J LL D L {E F )] 2 [J RR D R (E F )] 2 
2 2 
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= 271 



[J LL D L (E F )] 2 , [J RR D R (E F )] 2 



J 2 LR D L (E F )D R (E F ) 
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(J Z LL D L (E F )D R (E F ) 



dE f 
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- J 2 LR D L (E F )D L (E F )- 
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(16) 



The change in the density of states associated with the chem- 
ical potential gradient across the junction triggers a modifi- 
cation of the damping a that is linear in voltage. For typical 

L/R 

Fermi energy E F of the order of several electron-volts, the 
voltage dependence of a may become very notable. This volt- 
age driven effect may be expressed in terms of oco and oci (V) 
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with oc(V) = (Xo + 0Ci (Eq.(|2j). Here 

oo = n{j 2 LL [D(E L F )} 2 + 2J 2 m D(E L F )D{Ef) + J 2 RR [D(E R F )] 2 ) , 

1 / lj_ _Jo_\ 

Ul ~ 4e\T \ 2 [D(E^)+D(E§)}\[D(E^ + [D(E*)?) 

<WW!)(4^-4^) 

+ ^(P(^)]^-[ D (^)]»)](.17, 

CONCLUSIONS 

In conclusion, we present a theoretical study of Landau- 
Lifshitz-Gilbert damping (Eq.Q) for a localized spin S in a 
junction. The exchange interactions between the localized 
spin and tunneling electrons leads to additional dissipation of 
the spin motion, see Fig.Q. In the presence of an applied 
voltage bias V, the damping coefficient, i.e., Gilbert damping, 
is modified in linear order in V for the leads with particle-hole 
asymmetry in the Density of States. 

ACKNOWLEDGEMENTS 

Work at LANL was supported by the US DOE under LDRD 
X1WX. 

APPENDIX A: DERIVATION OF THE EFFECTIVE ACTION 

Here we will give a detailed derivation of the effective ac- 
tion for a spin. From Eq.©, we can extract a quadratic form 
of spins with the aid of the well known identity In detM 
Tr InM. In order to tabulate the expansion of Tr InM perturba- 
tively, we define matrices Mo and M\ , 



~ Ln=i n ' we can wr i te the effective action as 

Seff (S(x) ) ~ So + Tr InMo + Tr(M l Mi ) - ^Tr(M Q - L Mj f + ■ ■ ■ , 

(18) 

where So is the sum of the first and the second term in Eq.(|8jl. 
The third term in Eq.ljSJ (and consequent last term shown in 
Eq. (I18» is the first non-trivial contribution to the spin equa- 
tion of motion. Its evaluation is straightforward, 

Tt(Mq 1 M\) 2 

K 2 K\ 

(-/co' + ^ (5p )- 1 (7p a (co' - co)) / , 

where, repeated indices are implicitly summed over. Then, we 
find 

AS=^Iv(M Q 1 M 1 ) 2 

= 4 L L L L L + U)-\T a ^ Co') W 

P a>«>' a|5 oo' kea pep 

(-to' + ^-^Tpa^'-O)))^ 

= LLL LWpa^T^^-^-) -S(C0 m ) 

= -JdzJ c/t'S(t)-S(t / )2^(t-t')- 

Here, /f(x — x') denotes the integral kernel defined in 
Ea.lTrol, Upon invoking the identity 2S(x) -S(x') = ((S(x)) 2 + 
(S(x')) 2 ) - (S(x) -S(x')) 2 , the effective action becomes that 
ofEq.fl]]. 



M K p =(M )kp + (Mi)kp 

(M )kp = (( 'C0 + ^cd)5ap§oo')§KP APPENDIX B: THE DERIVATION OF THE SPECTRAL 

j DENSITY 

{Mi) K p =-^(r ap (co-fi)')W, 

We return to Ea. (llO> derived in Appendix A, and express 
where K = (k,(d) and P = (p,(a') with fermionic Matsubara the sum as a contour integral following standard procedures, 



frequencies co and co'. Employing the expansion ln(l +x) = e.g. 1 16], to obtain 
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K, 



ap 



keapefi 
'J n p,Ji 
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H[/fe)-/fe)]8(co+^-^ 



ieapep 

cosh(p/2_^]co 
aca/(a>) — . , ,„ — -— — 
o sinh(pco/2) 



P m — 1 1 — e-P m 

cosh(p/2- |t|)co 
sinh(pco/2) 



(19) 
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where /(©) denotes the spectral density in Caldeira-Leggett 
theory. The standard contour employed here is shown in 
Fig.Q. The symbol P in Eq.J19I> denotes the principal part 
of the integral. 




FIG. 3: The standard contour employed in Eq. <19> in order to evalu- 
ate the Matsubara sum of Eq.(lO|. The crosses along the imaginary 

axis denote the Matsubara frequencies. The spectra j densky ^ Eq ^ 



TZ poo r-<x> 

J{(a) = -f L d&ffo) /W^)[/fe)-/(^)]5(co + ^-^) 

** ^F F 

+ J -f f^d&Nfa) /^^(&)[ffo)-/(&)]8((D+k-&) 

+4 R r d&i<&) r j^ R N^ R )[f^ L )-f^ R m^+^ L -^ R ) 

~ = 0)N(% L = 0) + J -fN(^ = 0)N(% R = 0) + f iR N(k L = 0)N£ R = 0)) CD 

= (^fD L (E L F )D L (E L F ) + J -fD R (E R )D R (E R )+J 2 LR D L (E L F )D R (E R )) co, (20) 

I 

where D L I R (E F ^ R ) denotes the density of states at the Fermi If we apply a voltage leading to a chemical potential drop of 
energy level of the left/right lead. 
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E^-E§ = (E F + Afj L ) - (E F +A/j r )= eV, then Eq.O fol- 
lows. This, in turn, leads to Eq. (I16> . then Eq.( I13> follows. 



APPENDIX C: THE SPIN EQUATION OF MOTION 

If /(to) = then, from Eq. J12> . the non-local in time ker- 
nel of the action (Eq.(|5)) is K(t) ~ \. We thus obtain from 
Eq.O, 



(21) 



The functional derivative of ASdis with respect to S(x) is 



8A5 rffa a /• , (5(t)-5(tQ) 

&S(t) 7ti (x-x') 2 

From the free portion of the action (the first two terms of 
Eq.©), we have 



8So . 1 dS(z) -.. . r 



5S(x) S 2 £/x 



(23) 



Adding Eqs. (I22I23I >. equating the sum to zero, cross multiply- 
ing with S(x), and changing x — ► it, we obtain Eq.Q. 
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